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Abstract 

We construct strong solutions for a nonlinear wave equation for a thin vibrat- 
ing plate described by nonlinear elastodynamics. For sufficiently small thickness 
we obtain existence of strong solutions for large times under appropriate scaling 
of the initial values such that the limit system as h — > is either the nonlin- 
ear von Karman plate equation or the linear fourth order Germain-Lagrange 
equation. In the case of the linear Germain-Lagrange equation we even ob- 
tain a convergence rate of the three-dimensional solution to the solution of the 
two-dimensional linear plate equation. 
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1 Introduction 

In the present contribution we study the nonlinear wave equation for a thin vibrating 
plate (or beam if d = 2). The plate is assumed to be of small but positive thickness 
h > and satisfies the equations of three-dimensional nonlinear elastodynamics. 

In order to explain the result and the model under consideration, let us start by 
recalling some facts and results for the corresponding variational problems, see [7] 
for further details. We consider the elastic energy 

E h (z) = j-f (w(Vz{x)) - f h ■ (z(x) - x)) dx, 

where fi^ = Q' x (— \) is the reference configuration of the thin plate, Q! C M d_1 , 
d = 2, 3, is a suitable bounded domain, and z : —* M rf is the deformation of the 
plate. For simplicity, we will restrict ourselves to the case d = 3 in this introduction. 
Rescaling to O = O' x (— |, ^), we obtain the rescaled energy 



E h (y) = / W(V h y(x)) -f h -\ y(x) 



where y(x) = z(x',hxs) with x' = (xi,X2) and = (d Xl ,d X2 , j^d X3 ). The limit as 
h — > depends on the asymptotic behaviour of f h . More precisely, let / be of order 
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h a . If a = 2 , then the energy E h is of order h@ with (3 = 2. The rescaled energy 
jpE h converges as h — > to the elastic energy from the geometrically fully nonlinear 
Kirchhoff theory in the sense of T-convergence. To the authors' knowledge there are 
no results on existence of solutions for the corresponding dynamic wave equation 
or on regularity of non-minimizing equilibria. Indeed even the precise definition of 
equilibrium is not completely clear since the isometry constraint Vy T Vy = Id for 
the limit map y: Q' — > M 3 makes the problem very rigid; see Hornung [9, 10] for 
recent progress. If a > 2 and (3 = 2a — 2, then the limit energy can be described as 

where e(U) = sym (VJ7) , 

U = lim^ff^-M'), V = ^yl (1.1) 




2(q - 2) if 2 < a < 3 
a — 1 if a > 3 



(1.2) 



where Id'(s) = (x 1 ,x 2 ) T and Q 2 : ^ 2x2 -> R is related to Q 3 (F) := Z? 2 W(Id)(F, F) 
by 

Q2(G) = min Qs(G + a (g) e 3 + e 3 <g> a). 

aeR 3 

Here 

'+oo if2<a<3, 
A a = < 1 if a = 3, 
s ifa>3. 

Thus for 2 < a < 3 one has the "geometrically linear" constraint 2e{U)+W®W = 
0, which again has so far prevented the rigorous study of the associated dynamic 
wave equation or non- minimizing equilibria. For a = 3 (and therefore (3 = 4) one 
obtains the von Karman plate theory and for a > 3 (and therefore f3 > 4) one 
obtains a linear Euler-Lagrange equation (linear Germain-Lagrange theory), which 
for isotropic materials reduces to the biharmonic equation. 

Here we study the cases a = 3, [3 = 4 and a > 3, (3 = 2a — 2 > 4 in the dynamic 
situation. The equations of elastodynamics arise from the Lagrangian 

I j n {^f ~ W(Vz(x)) + A) dx = - W(V h y(x)) + /\) 

and solutions formally preserve the total energy 

(^f + W(V h y(x)) - f h y^j dx (1.3) 

In view of (1.1)-(1.2) we expect that 

y 3 ~ /i, - Id' ~ /j 2 for a = 3, /3 = 4 

y 3 ~ /i°- 2 , M - Id' ~ /i"- 1 for a > 3, /3 = 2a - 2 > 4 
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The idea to balance the kinetic and potential energy in (1.3) suggests to rescale time 
as r = ht if a = 3. Then the total energy becomes 

£ « = A4 I(^ + ^™ I), -Il)' b 

and with fh = h~ 3 f!^es the evolution equations is 

^d 2 T y - ±div h DW(V h y) = 

or equivalently 

d 2 T y - ^div h DW(V h y) = hf h , (1.4) 

where fh ~ 1 as h — > . In the case a = 3 we will show existence of strong solutions 
of (1.4) for well-prepared and small data in a natural scaling with respect to h and 
time r G (0,To) with To > sufficiently small. In particular we assume that the 
rescaled fh is small, cf. Section 3.1 below. - Note that the small time interval 
(0, To) for r turns over to a large time interval (0, TqH~ 1 ) in the original times scale 
for t. In the case a > 3, we will use the same time scale. Then we are able to show 
existence of strong solutions for r G (0, T) for any T > provided that fh ~ h a ~ 3 
and suitable initial data, cf. Section 3.1 below. We note that this time scale is 
subcritical for this scaling of fh- In this case we are even able to construct the 
leading term of the solution y = yh as h — > provided W(F) = dist(T, SO(3)) 2 , 
cf. Section 4. 

Together with [1] this shows that after the natural time rescaling and for well pre- 
pared data of the correct size solutions of the 3-d nonlinear elastodynamics converge 
to solutions of the dynamic von Karman equation or linear von Karman equation 
depending on the size of the data. We note that a similar result in the case of sta- 
tionary solutions was shown by Monneau [18] if the limit system are the von Karman 
plate equations. Ge, Kruse and Marsden [8] have taken an alternative and very gen- 
eral approach to study the limit from three-dimensional elasticity to shells and rods 
by establishing convergence of the underlying Hamiltionian structure. This suggests, 
but does not prove the convergence of the corresponding dynamical problems (see 
e.g. recent work by Mielke [17] for the question on the relation of the convergence 
of the Hamiltonian and the convergence of the resulting dynamical problems) . Gen- 
eral information and many further references on the dynamics of lower-dimensional 
nonlinear elastic structures can be found in the book by Antman [3]. For results on 
existence of weak and strong solutions of the non-stationary von Karman plate equa- 
tions we refer to e.g. Chen and Wahl [5], Koch and Lasiecka [12], Lasiecka [15], Koch 
and Stahel [13]. For a survey on results and open problem of nonlinear elasticity, 
stationary and non-stationary, we refer to Ball [4]. 

Let us explain the strategy of our proof and the main difficulties. Basically, the 
strong solutions are constructed by the usual energy method as presented e.g. in 
the books by Majda [16] and Dafermos [6]. (For a more abstract and general version 
see e.g. the classical paper by Hughes et al. [11].) The starting point in the method 
is the conservation of energy: 

dty(t)\\h { n) + ^J n w{ y y y) dx + j a h fh ■ y(t) dx ) = 
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which follows from (1.4) by multiplication with dty under appropriate boundary 
conditions. (Here and in the following we replace r by t.) Moreover, differentiating 
(1.4) with respect to x one gets a control of 



d 
dt 



^ll^y(*)lli 2( n) + p J a D 2 W(V h y)d%V h y : d%V h y dx) = R p , (1.5) 



where the remainder term Rp can be controlled with the aid of the Gronwall in- 



equality once the left hand side controls dxVhy suitably. To this end it is essential 
to have the coercive estimate 



1 



D 2 W{V h y)V h w : V ' h w dx > c 



h 



£h{w) 



L 2 (n) 



(1.6) 



where £h(w) = sym(V hw), cf. (3.34) below. By Korn's inequality in the present 
h -dependent version we have 



IIVhHI^n) ^ c 



cf. Lemma 2.1 below. Therefore we will have one order of h better decay of the 
symmetric part of V^y than for the full gradient/the skew-symmetric part. To 
obtain (1.6) (and similar estimates) it will be essential that 



^lkh(y) - /||l<» + II 



\L°° 2i 



< eh 



for some sufficiently small e > and to treat the symmetric and asymmetric part 
carefully in a Taylor expansion of D 2 W(Vhy) around /, cf. Sections 2 and 3.3 for 
the details. 

Several technical difficulties arise from the fact that we are dealing with natural 
boundary conditions at the upper and lower boundary x& = ±i. In tangential 
direction we assume periodic boundary conditions. First of all, in this situation it is 
easy to differentiate in tangential and temporal direction to obtain (1.5) with dxiv 
replaced by d^w, where z = (x', t) and x' = (x\, . . . , x^-i) ■ Therefore we are using 
anisotropic L 2 -Sobolev spaces of sufficiently high order to control V^y in L°° . In 
particular, one of the basic spaces is 

V(Q) = {«£ L 2 (Q) : Vu,d Xj Vu e L 2 (ft), j = l,...,d-l} ^ L°°(ft) 

if d = 2,3. Note that V(Q) is slightly larger than H 2 (9) and that u £ H 2 (tt) if and 
only if u € V(Q) and d 2 d u G L 2 (Q). Moreover, since we are dealing with natural 
boundary conditions, we want to keep the equation in divergence form. Therefore 
we do not use the identity 

div^I^y) = D 2 W(V h y) ■ V 2 y 

to obtain a quasi- linear system. Instead we differentiate (1.4) with respect to time 
and solve 

8 t w - ^div h (D 2 W(V h y)V h w) = hd t f, 
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where w = dty ■ Unfortunately we cannot solve the latter equation directly since 
we are missing suitable control of d x y , j = 1, . . . , d — 1. Therefore we first replace 
D 2 W(Vhy) by suitably smoothed coefficients A n (Vhy) and construct a solution 
w n , y n , respectively, for the smoothed systems. Once we have one solution y n at 
hand, we can differentiate it with respect to Xj , j = 1, . . . , d— 1 and get a solutions 
of 

8 t w] - -^dw h A n (V h y n )V h w™ = hd Xj f, 

where = d Xj y n . These solutions satisfy the same estimates as w n = dty n 
(uniformly in < h < ho and the smoothing parameter n G N). Then we can pass 
to the limit n — > oo to obtain a solution of the original system. 

The structure of the article is as follows: In Section 2 we introduce some notation 
and derive some preliminary results. Our main result is presented in Section 3.1. 
Afterwards we introduce our approximate system used for the construction of strong 
solutions in Section 3.2. The essential results for the linearized system are derived in 
Section 3.3, which are applied in Section 3.4 to obtain a strong solution first locally 
in time for fixed n G N. Then uniform bounds in J 1 , ft and n G N are shown in 
Section 3.5 and our main result is shown by first extending the solution for small 
times in [0, T] for some T > given and then passing to the limit n — > oo . Finally, 
in Section 4 we derive a first order asymptotic expansion as h — > in the case that 
the limit system is linear, i.e., (3 > 4, and W(F) = dist(F, SO(d)) 2 . 
Acknowledgements: This work was partially supported by GNAMPA, through 
the project "Problemi di riduzione di dimensione per strutture elastiche sottili" 2008. 

2 Notation and Preliminaries 

For any measurable set M C M> N the inner product of L 2 (M) (w.r.t. to Lebesgue 
measure) is denoted by (.,.)m- Moreover, H k (Q), k G No, denotes the usual I? - 
Sobolev spaces. If A is a Banach space, then the vector- valued variants of L 2 (M) 
and H k (M) are denoted by L 2 (M;X), H k (M; X) , respectively. Furthermore, 
C fc ([0,T];A), k G N , denotes the space of all k -times continuously differentiable 
functions /: [0, T] — ► X . 

For the following Q = (— L, L) d_1 x {-\,\), St = {—L,L) d ~ 1 , d = 2,3, x = 
{x',x d ), where x' G R*" 1 , let V h = (V^, \d Xd ) T , V x , t = (d t ,V x ) and let 

e h (w) = sym(V^), e(w)=e 1 (w), 

if w: M C R d — > M. d is a suitable vector field. Here symA = \{A + A T ) and we 
denote skew A := \{A — A T ). Moreover, we denote z = (t,x'), where zq = t and 
Zj = Xj for j = 1, . . . , d — 1 . 

For s > 0, s G" No, we define L 2 -Bessel potential spaces 

H S (Q) = {/ G L 2 (Q) : f = F\ n for some F G H s (R d )} 

as usual by restriction, equipped with the quotient norm. Since £1 is a Lipschitz 
domain, there is a continuous extension operator E such that E: H k (VL) — ► H k (R d ) 
for all k G N, cf. Stein [20, Chapter VI, Section 3.2]. Hence H S (Q), s > 0, is 
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retract of H s (M. d ) and we obtain the usual interpolation properties, cf. e.g. [21]. In 
particular, we have 

(H s »(n),H s ^n))e,2 = H s (n), a = (1 -o)s + o Sl , (2.1) 

for all 9 G (0, 1) , s > 0, where (., .)q <p denotes the real interpolation method. 

If < T < oo and X is a Banach space, then BUC([0, T]; X) is the space of 
all bounded and uniformly continuous functions /: [0, T) — > X . Now let Xo,Xi be 
Banach spaces such that X\ Ao densely. Then 

^(0, T; X ) n L"(0, T; X l ) - BUC([0, T]; (X ,X 1 ) 1 _^ p ) (2.2) 

for all 1 < p < oo continuously, cf. Amann [2, Chapter III, Theorem 4.10.2]. If 
Xq = H is a Hilbert space and H is identified with its dual, then X\ ^ H ^ X[ 
and 

\j t \\ft H = (j t f(t)J(t)) x , tXl for almost all t G [0,T] (2.3) 

provided that / G L p (0,T;X!) and fj G L p ' (0,T; X[) , 1< p < oo, cf. Zeidler [23, 
Proposition 23.23]. In particular, (2.3) implies 

sup \\f(t)\\ 2 H < 2 (||aj|| L2(0iT;X;) ||/|| L2(0iT;Xl) + ||/(0)||^) . (2.4) 

te[o,T] v ' 7 

Replacing f(t) by t/(t) and (T — t)f(T — t), one easily derives from the latter 
estimate 

sup \\f(t)\\ H < Qrll/llia^jjll/lli,^) (2-5) 
*e[o,T] v > - iy v > > v 

for some Ct > depending on T > . 

In the following £ n (y), n G N, denotes the space of all n-linear mappings 
A : V n —>■ M for a vector space V. Moreover, if ^4 G £ n (V) , n > 2, and xi, . . . , G 
V , 1 < k < n, then ^4[xi, . . . , x^] G C n ~ k {V) is defined by A[x\, . . . , Xfc](xjt+i, . . . , x n ) = 
A(xi, ... ,x n ) for all x k+1 , . . . ,x n G V . 

We introduce the scaled inner product 

A : h B = -^sym^l : symB + skewyl : skewS, A,B G M. dxd ,0 < h < 1, 

and lAj/j = vG? :h A where A : B = J2i j=i a ijbij ■ This choice of inner product is 
motivated by the Korn inequality in thin domains, see Lemma 2.1 below. Of course, 
:i coincides with the usual inner product : on 'R dxd and therefore \A\\ = \A\. For 
W G £ n (R dxd ) we define the induced scaled norm by 

\W\ h = sup \W(A 1: ...,A n )\. 

\Aj\ h <l,j=l,...,n 

Note that, since \A\ h > |A|i = \A\ for all A G M. dxd , we have \W\ h < \W\i =: \W\ 
for any W G C n (R dxd ) and < h < 1. 

As usual we identify £ 1 (M dxd ) = (M dxd )' with ]R a!xd . But one has to be careful 
whether this representation is taken with respect to the usual scalar product : on 
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R dxd Qr with respect to i e ^ G £i( R dxd) is identified with ^ G R dxd such 
that 

= A : h B for all B G M dxd . 

If nothing else is mentioned, we identify (R dxd )' and R dxd using the standard 
inner product :. In particular, if W G C 1 (i7), f7 C M dxd and A G U, then 
£>W(A) G (R dxd )' 9* lR fl!xd coincides with 



DW(A) : B = jW{A + tB) 



for all B G M dx,i . 



t=o 

Furthermore, VF G £ 2 (R dxd ) is usually identified with the linear mapping IR rfx<i 
R dxd d e g ne d ^ 

WA : B = W(A, B) for all A,B £ R dxd . 
Finally, we denote by 

WWLliM-c^^i)) = \\ w \\lI(m) = {^J m \W{x)\ p h dx^ 

if 1 < p < oo and with the obvious modifications if p = oo. Here M C R' 
is measurable. Moreover, for / G L p (M;R dxd ) the scaled norm ||/||£P(M;R dxd ) = 
\l p (m) 1S defined in the same way. 

We now state the relevant Korn inequality in thin domains. 



Lemma 2.1 There is a constant C such that 



l|Vh«|U 2( n) < C 



\e h {u) 



(2.6) 

L 2 (n) 



for all < h < 1 and u G H l (Q) d such that u\ Xj =-L = u\x j =L , J = l,...,d — 1- 

Proof: For clamped boundary conditions the Korn inequality in thin domains was 
proved by Kohn and Vogelius [14, Prop. 4.1]. They mention that the result also 
holds without boundary conditions, modulo infinitesimal rigid motions. For the 
convenience of the reader we provide a proof of Lemma 2.1. 

First we prove the case d = 2. Let := (— L, L) d ~ l x (— |, |) and let u G 
H l (Q,h; R 2 ) satisfy the boundary conditions u\ Xj= -L = u\ Xj =L, j = 1, ...,d— 1. 
First of all by a simple scaling in x^, (2.6) is equivalent to 

C 

\WA^(Q h ) < ■^■||(V«) awm ||i ( 2(n h ) (2-7) 

Let N h be the integer part of ^ and let 1^ := j^- . We set := {— L + : 
k = 0, . . . , Nh — 1}. By applying Korn inequality on the set (a, a + £h) x (— §, |) for 
every a 6 4, we can construct a piecewise constant function ^4: (— L,L) — ► M 2x2 
such that A(xd) is skew-symmetric and 

f \Vu-A\ 2 dx<C [ \e{u)\ 2 dx. (2.8) 
Jn h Jn h 
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Note that, since ^ is bounded from above and from below, we can use the same 
Korn inequality constant on each set (a, a + 4) x (— |, |) . 
We claim that 



/ \A{x\) - A \ 2 dx < ^ f \e(u)\ 2 dx. 
Jn h h J Uh 



(2.9) 



where A := A(-L). 

Let us fix a G Jh and let 6 := a + A4> with A G {0,1}. By applying Korn 
inequality on the set (a, a + 24) x (— |, |) we have that there exists A G M 2x2 such 
that 

|e(n)| 2 dx. 



Ox(-f.f) 



/ \Vu-A\dx<C [ 

J (a,a+2e h )x (-|,|) J(a,a+2£ h ) 

From this inequality we deduce 

/i4|^(fr)-i| 2 < 2/ |Vu- A(xi)| 2 dx 

y( 6 , 6+ 4)x(-f,f) 

>/ 

7(6 

■/ 



+2 / |Vu-i| 2 dx 
/(6,6+4)x (-§,§) 



< c 



'(a,a+24)x(-|,|) 



|e(«)| 2 dx. 



Combining the previous inequality for A = and A = 1 , we obtain 
M h \A(a) - A(b)\ 2 < 2M h (\A(a)-A\ 2 + \A(b)-A\ 2 ) 



< C 



(a,a+24)x 



\e(u)\ 2 dx. 



As A is constant on each interval (a, a + 4) , this is equivalent to say that 

/ \A( Xl +4) -A{ Xl )\ 2 dx < C [ \e(u)\ 2 dx. (2.10) 

7(a,a+£ h )x(-f,f) J(a,a+2^)x(-|,f) 

Let us set Ikj ■= —L + £h(k, k+j) . By convexity we have the following estimate: 

N h -1 

/ \A(x 1 ) - A \ 2 dx = h ^ / \M x i) ~ A o\ 2 dxi 

N h -1 k~l 

= h Yl / I H {A{x 1 -mi h )-A{x 1 -{m + l)i h j) 

k=0 Jfc .! m=0 
N h -1 k-1 

By (2.10) we deduce 

iV/j — 1 fc-l 



<ixi . 



e(n)| dx. 
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It is easy to see that for every k = 0, . . . , — 1 

k-l 



E 



m=0 U 1 k-m-l,2 x (~^>f) 



e(u)| 2 dx < 2 / |e(u)| 2 dx. 



Therefore, we conclude that 

! \A( Xl )-A \ 2 dx<CNlf \e(u)\ 2 dx, 

which proves claim (2.9). 

Combining (2.8) and (2.9), we conclude that for every u G H 1 ^^,^ 2 ) there 
exists a constant skew-symmetric Aq £ M 2x2 such that 



\Vu - A \ 2 dx < 



C_ 



\e(u)\ 2 dx. 



Since 



[ jf—- [ (skwVu) dx - A dx < [ |(skwVn) - A \ 2 dx, 
Jn h \^h\ Jn h Ju h 



we also have that 



Vu 



\^h\ Jn h 



(skwVu) 



dx < 



C_ 



f Hut 



dx 



(2.11) 



for every u G H 1 ^^, M 2 ) . 

Now, if u is periodic in tangential direction, then 



/ TFT7 / ( skwVu ) dx = [ / d 2 ui 

Jn h Jn h Jn h Jn h 

1 



dx 



\n 



h\ Jn h 



(d 2 ui + diu 2 ) 



dx 



< I Hut 



dx, 



which, together with (2.11), provides us with the desired inequality. 
In order to prove the case d = 3, we use that (2.6) for d = 2 implies 



0, 



l 

h" x 3 



d X:} J \U 3 



< 



C 



L 2 (n) 



\dx 3 ) \u 3 



sym 



c 



L 2 (Q) 



< -^\\(^hU) sym \\ L 2 (n) 



for j = 1,2 and any u £ H (CI) . Moreover, applying Korn's inequality in (— L, L)' 
with periodic boundary conditions, we obtain 

W^x'u'Wl'I^ < C||(V x'u') S ym\\L 2 (n) < C\\(\7 x U) sym \\ L 2(ty, 

where v! = [u\,U2) T ■ Altogether this proves (2.6) for d = 3. ■ 
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Remark 2.2 The latter lemma shows that \\^h(u)\\ L 2fQ\ is equivalent to \\Vhu\\ L 2^ 
with constants independent of < h < 1 . 

We denote by H^Q) the space Hp er (Q) d n {u : J n u(x) dx = 0} equipped with the 

L 2 (f2) 



norm 



h 



u G ff 1 ^ 



and H h its dual space with norm 

\H-\n) = SU P | </. v) h -\hI '■ u G with 



L 2 



Furthermore, we denote 

H£r(n) = {/ G # m (tt) : ^/U, = _ L = d2f\ Xj=L ,j = 1, . . . ,d - 1, |a| < m - 1} . 

Throughout this contribution the following anisotropic variant of H™ r (Q) will 
be important: 

H mi > m *(n) = { U £L 2 (n):V k x/ d l Xd ueL 2 (n),k = 0,..., mi ,l = 0,...,m 2 , 

dx'd l Xd u\ Xj =- L = d^d l Xd u\ Xj=L ,j < d- l,\a\ < mi - 1,1 < m 2 } 
where mi G N, m 2 G No . The spaces are equipped with the inner product 

|a|<rrti,fc=0,...,TO2 

Please note that periodic boundary conditions are included in the spaces H mi ' m2 (tt) 
in contrast to the space H m {£l) , where we denote them by a subscript "per" in order 
to be consistent with the usual definition of H m (Q) . 
Similarly, an anisotropic variant of L p will be useful: 

L™(fl) = {uin-vRiHari.Oll^.i^eLPa-^L)"- 1 )} 

where 1 < p, q < oo equipped with the norm 

||/|| ^ = ll ||n(xi '- )|| ^(-ii)IL,((- w - 1) - 

We note that from the usual Holder inequality it follows that 

ll/s , llLp.«(n) < ||/||LPi.9i(n)ll5llLP2.92(n)) 
for all 1 < pi, qi,P2, q2 < oo such that 



1 1 1 

- = — + —: 

P Pi P2 



1 1 1 
- = — + —. 
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Lemma 2.3 Let d = 2,3. Then 

H 1 ' ^) ^ L p ' 2 (n), H 2 '°(Q) ^ L°°' 2 (ft), H 1 ^) ^ L 4 '°°(ft) 
continuously for p = oo if d = 2 and any 1 < p < oo if d = 3 . Finally, let 

V(Q) := H 1 ' 1 (Q)nH 2 '°{n). 
Then V(Q) (7° (17) continuously. 

Proof: The first embedding follows from H 1 ^') <^-> L p (f2') and the second from 
# 2 (ft') ^ L°°(0') since d = 2,3 and fi' = (— L,L) d_1 . The third embedding 
follows from 

^(-i, 1; L 2 (tt')) n L 2 (-I, I; ff 1 ^')) - Bl7C([-£, if* (O')) 

and H^{Q!) <^-> L 4 (lV). Finally, the last embedding follows from 

L 2 ( _ i > i . ff i+*((_ L> L)^ 1 )) n ^(-I, I; ^((-L, L)^ 1 )) 
^ BUC([-l \];H^ ((—L, Lf- 1 )) - C°(fi) 

where k = d — 2 because of (2.2) and Sobolev embeddings. ■ 

Remark 2.4 The spaces H 1,0 (Q) and V(f2) are the fundamental spaces, which will 
be used to solve the evolution equation. We note that 

fev(n) & f,vfeH^°(n). 

Most of the time we will estimate f G V(Q) by the h -dependent norm 

\\f\\v h := \\(f,V h f)\\ m ,o {n) . 

Because of the embedding V(Q) <^-> L°°(Q), we are able to show that V(Q) is an 
algebra with respect to point-wise multiplication. More precisely, we obtain: 

Corollary 2.5 Let d = 2,3. Then there is some C = C(fi) > such that 

|| (mi • u 2 ,V h (u 1 ■ u 2 )||tfi,o(n) < C\\(ui,V h ui)\\ H i,o( n )\\(u 2 ,V h u 2 )\\mfi(n) ( 2 -12) 

for all u±,u 2 G V(tt) uniformly in < h < 1. Moreover, if F G C 2 (U) for some 
open U C l w , iV G N, and u G V^fi)^, i/ten /or ewery i? > i/tere zs some C(i?) 
independent of u such that 

\\(F(u),V h F(u))\\ Hl ,o {n) <C(R) if\\(u,V h u)\\ H i,o {n) <R (2.13) 

uniformly in < h < 1 and if u(x) G C7 /or all x G f2 . 
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Proof: First of all (2.12) follows from (2.13) by first considering ||«i||y h , ||«2||y h < 1 
and F(ui,u 2 ) = u± ■ u 2 together with a scaling argument. 
Hence it only remains to prove (2.13). First of all, 



d Xj F(u) 



d Xj d Xk F(u) 



DF{u)d Xj u 

DF{u)d X] d Xk u + D 2 F(u)(d Xj u, d Xk u) 



where DF(u), D 2 F(u) are uniformly bounded since u G C°(fi) and u(x) G U for all 
x G S7. Therefore VhF(u) G L 2 (17) can be easily estimated. Hence it only remains 
to consider the second order derivatives. To this end we use that 



1 



D F(u)(d Xj u, jS Xd u) 



< C\\d Xj u\\ L 4,<=o(p-) 



L 2 (Q) n L 4 > 2 (Q) 

<C'(R)\\(u,V h u)\\ m , om 



-rd Xj u 
h d 



< C\\d Xj u\\ H i(ci) 

Hi.0(O) 

for all j = 1, . . . , d — 1 due to Lemma 2.3. Similarly, 

\\D 2 F(u)(d Xj u,d Xk u)\\ LHn) < C'(R)\\(u,V h u)\\ m , 0{Q) 

for all j, k = 1, . . . , d — 1. From these estimates the statement of the corollary easily 
follows. ■ 

For the following let W: B r (I) C M. dxd — > R be a smooth function for some 
r > which is frame invariant, i.e., W(RF) = W(F) for every F G M dxd — > IR dxd 
and i? G 50(d) , and such that DW(I) = and D^W(I) : M dxd -► M rfx<i is positive 
definite on symmetric matrices. Moreover, we set W{G) = W{I+G) . The estimates 
of derivatives of D 2 W(Vhu) will be essential for the proof of our main result and 
will be based on the following lemma: 



Lemma 2.6 There is some constant C > 0, e > 0, and A G C°°{B e (0); C 3 (R dxd )) 
such that for all G G R dxd with \G\ < e we have 



D 3 W(G) = D 3 W(0) + A(G), 



where 



\D 3 W{0)\ h < Ch for allO < h < 1, 
\A{G)\<C\G\ forall\G\<e. 

Proof: First of all, if |G| < e for e > sufficiently small, we can use a polar 
decomposition I + G = RU , where R G SO(d) and U is symmetric and positive 
definite such that U 2 = (I + G) T (I + G) . From frame invariance we conclude 
that W(I + G) = W(U) = W(U 2 ) = W(I + 2symG + G T G) for some smooth 
W: V C M. dxd — ► R, where V is some open neighborhood of /. For this proof we 
denote A s = sym A . Straight-forward calculations yield 

DW{F){H) = DW(U 2 )(2H S + H T G + G T H) 
D 2 W(F)(Hi, H 2 ) = D 2 W(U 2 )(2H 1<s + HfG + G T H 1 ,2H 2 , s + HjG + G T H 2 ) 

+DW(U 2 )(HfH 2 + H^Hr) 
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and 



D 3 W(F)(H 1 ,H 2 ,H 3 ) = 

D 3 W{U 2 ){2H lyS + HfG + G T Hi, 2H 2 , S + H%G + G T H 2 , 2H 3>S + H^G + G T H 3 
+D 2 W(U 2 )(HlH 2 + H 2 Hi, 2H 3 , S + H^G + G T H 3 ) 
+D 2 W(U 2 )(HlH 3 + flffTi, 2F 2jS + H 2 G + G T H 2 ) 
+D 2 W(U 2 )(H 2 T H 3 + #J # 2 , 2# M + flf G + G T Hi) 



where F = I + G . From the latter identities the statements immediately follow. 
For the following we denote 



■■= ( E 11^^11^(01 

\a\<mi,j=0,...,m 2 

l 

\ A \\h? ■= ( E H^lli^n) 

v |a|<m 

where m,m 1 ,m 2 G N and 4 G fP ni > ma (fi) dxd , 4 G H m (n) dxd , respectively. 
Corollary 2.7 There are some e, C > smc/i £/iai 

p^^Xy^yOlUi^ 

Hjj' ) ||^2||x,2(n) ll^3|li,2(n) (2-14) 
for all y G y(0) dXfl! ,y2,y3 G L 2 (Q) dxd , 0<h<l and \\Z\\ La o {n) < mm(e,h) and 
\\D 3 W{Z){Y u Y 2 ,Y 3 )\\ L , {n) 

< ch i|Yiii H i (n) i|y2ii H i,o (n) i|y 3 || L 2 (Q) (2.i5) 

/or a// y G Jf 1 (n) dxd ,y 2 G H 1 >°(tt) dxd , Y 3 G L 2 (n) dxd , < /i < 1 and Z G 
L°°(ft) dxd with \\Z\\ L oo {n) < mm{e,h). 

Proof: The statement follows directly from Lemma 2.6, Korn's inequality due to 
Lemma 2.1, and Lemma 2.3. ■ 



3 Long-Time Existence for Thin Sticks/Plates 

3.1 Main Result 

We consider 

8 2 u h - ^dw h DW(V h u h ) = f h h 1+e in ft x / (3.1) 

13 



where W(G) = W{I+G), U = (— L, x (— |, |) , /3 = 4+20, which is equivalent 

to = a — 3, and I = [0, T*] for some T* together with the initial and boundary 
conditions 



DW(V h u h )e d 



Uh\ Xj=L 
(u h ,d t u h )\ t=0 







u h\ Xj =-L> 
(uo,h,ui,h)- 



l,...,d-l, 



(3.2) 

(3.3) 
(3.4) 



Here we assume that W : -B r (I) — > M is a smooth function for some r > which is 
frame invariant, i.e., W(RF) = W(F) for every F G R dxd -► M dxd and G SO(d), 
and such that DW(I) = and D 2 W(I): R dxd -► R dX(i is positive definite on 
symmetric matrices. - Note that the latter condition implies that D 2 W(I) is elliptic 
in the sense of Legendre-Hadamard: 



(D 2 W(I)a <g> 6) : a <g> 6 > c |a| 2 |6| 2 for all a,k 



(3.5) 



for some Co > 0. 

Theorem 3.1 Let 9 > 0, < T < oo, let f h G WftO, T; if 1 - ) n 14^/(0, T; i/ 2 >°) , 

< ft, < 1, and let u , h G H l > 3 (n) d n H 4 >°(n) d , u hh G H l > 2 (n) d n H 3 > (ti) d such 
that 



and 



DW{V h u , h )e d 



x d =±- 



D 2 W(V h u 0ih )V h u 1:h e d 



max 
fe=0,i 



^£/i(//i|t=o)>3t//i|t=0 



-J^h{ui+k,h),U2+k,h 



H 2 - k '°{VL) 



< Mh 1+e (3.6) 



+ max||^/ h || w i (0TH i, 0) < M (3.7) 



uniformly in < /i < 1 , w/iere 



/* 1+e M(=o + ^div ft DW(V fc «o,h), 



(3.8) 
(3.9) 



If 9 > 0, then there is some ho G (0, 1] and C depending on M such that for every 
0<h<h there is a unique solution u h G C 3 ([0, T]; iJ 1 ' ) n C°([0, T]; i/ 2 > 2 n i^ 4 ' ) 
o/ (S.1)-(S.4) satisfying 



max 

M<1 



d?d2u h ,V x>t d]-e h (u h ) 



Ch l+e 



C([0,T];ifi.°) 



(3.10) 



uniformly in < h < ho- If 9 = 0, i/te same statement holds with ho = 1 provided 
that M, T > ore sufficiently small. 
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3.2 Approximate System 

In order to construct a solution to (3.1)-(3.4), we first construct solutions for an 
approximate system. 

To this end, let P% : L 2 (Q) — ► L 2 (£l) , n € N, be the projection on the eigenspaces 
of the eigenvalues not exceeding n of — A^ with domain T>{— A^) = {u G H 2 er (Q) : 
d Xd u\ Xd=± i = 0}. Then P^ are orthogonal projections such that 

l|VfcP n Vlb(n) = ||(-A^)^P n /|| i2(n) = ||P„(-A ft )3/|| La(n) 

< ||(-A^)V|| L2(n) = ||V^/|| L2(Q) . (3.11) 

The same is also true for H m '°(Q) since — A^ and therefore also P% and (— A^)^ 
commute with tangential derivative d^, . Furthermore, the previous estimates imply 

\\d Xd PHf \\mn) < /»||V h /|| La( n) < ||V/|| L2( n) 

and therefore 

l|Vi^/|| La( n)<2||V/|| La(n) . 

Moreover, by standard elliptic theory all eigenfunctions are smooth. Therefore 
Tl{P%) C C°°(n). Finally, there is a constant C n > such that 

ll^/llfl*(n) < C n \\f\\ LHn) (3.12) 

uniformly in < h < 1 since 

\\PZf\\H*(n) < c(||/|| L2(Q) + ||AP n V|| 

+ ||A^/|| L2 (Q)) 



C 



< 



C (ll/H^n) + n\\PHf\\ L2(n) ) < C n \\f\\ LHQ) 



where we have used that ||A^ej|| i2 (Q) = Aj||ej|| i2 (Q) < ra||ej||£2(n) f° r each eigen- 
function ej to some eigenvalue Xj < n and 

Ci||A^|| L2(n) < HAaj/ullia^) + ^ IfeHll^O) - ° 2 W A HIl 2 (Q) 

for any < h < 1 . For notational simplicity we will write P„ instead of P„ in the 
following. 

To motivate the approximation, we note that 

DW(V h u) = D 2 W(0)V h u + [ D 3 W(TV h u)[V h u,V h u]{l-T)dT, 

Jo 

= D 2 W{0)V h u + F'{V h u), (3.13) 



by Taylor's expansion since DW(0) = DW(I) = 0. Therefore we define the approx- 
imations 

F n {V h u) := D 2 W(0)V h u + [ P n D 3 W{P n rV h u)[P n V h u, P n V h u](l - r) dr 

Jo 

= D 2 W(0)V h u + F^P n V h u) (3.14) 
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where n G N. Replacing DW(Vhu) by F„(V/ l u) in (3.1), we obtain the following 
approximate system. 

8 2 t u n h - ^dw h F n (V h v%) = f h h l+e infix (0, T) (3.15) 

F n (V h K)e d \ Xd=±1 _ = (3.16) 

u n h \ X]=L = u n h \ Xj= _ L , j = l,...,d-l, (3.17) 

K,«)l t=0 = K^<J- (3-18) 
Here UQ h ,u^ h will be chosen as solutions of 

«2,A = /t W A|(=o + ^div h F n (y h v% jh ) (3.19) 

«3,h = h 1+e d t f h \ t=0 + ^div h (DF n (V h ul h )V h ul h ) (3.20) 

together with J fi Uj h dx = j n Uj t h dx, j = 0,1, and the boundary conditions 

F n {V h ul h )e d \ Xd=±h = DF n (V h ul h )V h ul h e d \ Xd=±h = 0, (3.21) 

where ui^u-i^h are as in Theorem 3.1. We will show that (3.19)-(3.21) has a unique 
solution UQ h ,u™ h satisfying 



max 
fe=0,l 



^hiul th ),V^e h {ul h ),V 2 h ul h 



< C Mh 1+e (3.22) 

H 2 - k <°(n) 



for all < h < ho provided that ho £ (0, 1] is sufficiently small in the case 9 > 
and provided that M > is sufficiently small (and ho = 1) if 6 = 0. Here Co is 
some universal constant. 

Remark 3.2 If k = 0, the statement follows from Proposition 3.8 below. If k = 1 
and uo,h already constructed, the statement follows from Lemma 3.5 below and the 
Lemma of Lax-Milgram. 

The main step now consists in solving (3.15)-(3.18) under the same assumptions 
as in Theorem 3.1 and showing uniform bounds in n G N and < h < ho. More 
precisely, we show 

Theorem 3.3 Let 9 > 0, < T < oo , and let u 2t h, u 3,h, fh be as in Theorem 3.1 
. If 9 > 0, then there are ho £ (0,1] and C depending only on M,T such that 
for every < h < h and n G N there are unique solutions Uq w and u\ n of 
(3.19)-(3.21) and a unique solution u\ of (3.15)-(3.18) satisfying 



max 

M<1 



d 2 d2ul,V x , t d2\e h {ui: 



< Ch 1+e (3.23) 

C([0,T];//M>) 



uniformly in < h < ho , n G N . If 9 = , the same is true provided M, T > are 
sufficiently small. 
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Once Theorem 3.3 is proved, the main Theorem 3.1 is easily proved by passing to 
the limit n — > oo for a suitable subsequence using the uniform bounds due to (3.23). 

For simplicity we will write u n instead of u 1 ^ in the following. 

In order to construct a solution to (3.15)-(3.18), we differentiate (3.15)-(3.18) 
once with respect to time. This yields the system 

d 2 t w n - ^div h (A n (V h u n )V h w n ) = d t f h h 1+e infix (0, T), (3.24) 
A n (VhU n )V h w n e d \ Xd=± i = (3.25) 
w n \ Xj=L = w n \ Xj= _ L , j = l,...,d- 1, (3.26) 
(w n ,d t w n )\ t=0 = (w , n ,wi,n)- (3.27) 
for w n = dtu n = dtu\ , where 
A n (V h u n )V h w 

= D 2 W(0)V h w + [ P n D 3 W(P n TV h u n )[P n V h u n ,P n V h w]dr, (3.28) 
J o 

= D 2 W(0)V h w + P n (D 2 W{P n V h u n ) - D 2 W(0)^j P n V h w 

since F n (V h u n ) - D 2 W(0)V h u n = P n DW(P n V h u n ) - P n D 2 W(0)P n V h u n due to 
(3.13)-(3.14). - We note that A n (VhU n ) defines a symmetric operator on L 2 (Q) dxd 
since P* = P n and that 

DA n {V h u n )[V h v,V h w] = P n D 3 W(P n V h u n )[P n V h v,P n V h w}. 

Moreover, we have u>o,n = ui,n = u i h an< ^ 

wi,n = v% th = U2, h = fhh 1+e \t=o + ^div h (DW(V h u , h )), 

w 2 , n = ul h = u 3A = d t f h h 1+e \ t=0 + ^dw h (D 2 W(V h u 0)h )V h u hh ), 

provided (3.19)-(3.20) hold. Hence (wi n ,W2 n ) = (^1,^2) are independent of n G 
No and «o,nj^i,n- First we will solve (3.24)-(3.27) for w n and small times (depending 
on n G N) provided that < h < ho is sufficiently small ho € (0, 1] if 6 > and 
that M > is sufficiently small if 6 = (independent of n G N). Here u n is 
determined by w n via 

u n (x,t) = UQ h (x) + / w n (x,T)dT. (3.29) 
■/ 

Afterwards we will derive uniform bounds inO</i</io,nGN, and t G (0, T). 

Finally, we note that, if w n solves (3.24)-(3.27) and u n is determined by (3.29), 
then u n solves (3.15)-(3.18) since (3.24) implies 

d 2 u n - ^div h F n {V h u n ) = f h h 1+e + c, 

where c = c(x) is independent of t, and the initial condition dtw n \t=o = d 2 u n \ t =o = 
wi >n implies c = by the choice of wi >n . 
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3.3 Estimates for the Linearized Operator 



Recall that z = (t, x') with the convention that zq = t and Zj = Xj for j = 
1, . . . , d — 1 . Moreover, recall that V 2 = V^a/ = {dt, V x /) . Furthermore, P n = P% , 
n G N, < h < 1, denotes the smoothing operator defined above and we set 



Poo = I- 



Let Uh for some < h < 1 be given such that 



max 

|7|<1 



-^e h (dn z u h ),V x ^e h {Sn z u h ) 



Rh 



(3.30) 



where R G (0, i?o] f° r some < i?o < 1 to be determined later. We note that (3.30) 
implies that \e h (dlu h ) £ C([0,T]; V(fi)) and d?±e h (u h ) G C([0, T]; fl" 1 - ) . For the 
following we denote 

ll/llv h = ||(/,V ft /)|| H i,o, 
where / G V(O). Of course < \\f\\ Vh 



for all < h < 1 . 

Because of || V^Pn/H^ < ||V^/||^2, Korn's inequality (2.6), and since P ra com- 
mutes with derivatives with respect to z = (t, x') , (3.30) implies 



max 

M<1 



+ 



d]P n V h u h , d]P n -e h (u h ) 



d?P n V h u h , dfP n -e h (u h ) 



L°°(0,T;V) 



< dRh 



(3.31) 



L°°{Q,T;H 1 ' ) 

for some C\ > 1 depending only on the constant in the Korn inequality. 



Remark 3.4 TTte analysis in the following will be mainly based on (3.31 ). Therefore 
we will assume throughout this section that (3.31) holds for some given VhUh, n G 
No U {oo}, and < h < 1. - Of course, if Vh u h satisfies the stronger estimate 
(3.30), we will have (3.31 ) for any n G N U {oo} . 



Because of V(f2) L°°(fi), cf. Lemma 2.3, (3.31) implies in particular 



Pn^hUh,Pnj^e h (u h ) 



< MRh, 



L°°(0,T;L°°nV) 



(3.32) 



where M depends only on Q. Recall that W(A) = W{I + A) for all A € R dxd . 
In order to evaluate DW(P n Vf l U} l ) , we will assume that R$ > is so small that 
W G C°°(-BmRo(0)) and Mi?o < e, where e > is as in Corollary 2.7. 
Using (3.32) and (2.14), we obtain 



J (D 3 W(TP n V h u h (t))[P n V h u h (t),P n V h v},P n V h w 



< Co 



/I 



P n V h u h ,P n -e h (u h ) 



< C R 



~-£h(v) 



h 



V(P) 

L 2 (n) 



-£h(v) 



L 2 (n) 



L 2 (H) 

1 
h 



dr 



L 2 (n) 



(3.33) 
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uniformly in v,w G i^ er (ft) d , <t <T, < h <l. 
In particular, we derive 

^(A n (V h u h (t))V h v,V h v) L 2 {n) = ^(D 2 W(0)V h v,V h v) L 2 (n) 



+ 



J (D 3 W(TP n \7 h u h (t))[P n V h u h (t),P n \7 h v],P n V h v) L2 ^ dr 



> (c -CoRo) 



h £h{v) 



L 2 {Q) 



uniformly in v G H^ er (n) d , t G [0,T], 0<T<oo < R < R , < h < 1, where 
Co > depends only on D 2 VF(0) and f2. Hence, if Rq G (0,1] is sufficiently small, 
we have 



-^(A n (X7 h u h {t))X7 h v,X7 h v) L 2 {n) > ^ 



L 2 (Q) 



(3.34) 



for all v G Hl er {Vt) d , £ € [0, T] , < ft < 1, < R < R , and u h satisfying (3.31), 
where cq is as above and depends only on D 2 W(0) and Q. By the same kind of 
expansion for D 2 W and estimates one shows 



-j^(d Zo A n (y h u h (t))V h v,V h w) L 2 {n) 













< C'R 


\e h {v) 


L 2 




L 2 



(3.35) 



for all v,w G H^ er (Q) d , j = 0,...,d-l uniformly in < h < 1, t G [0,T], 
0<R<R , < T < oo. Therefore 



^2 (-4n(V^(t))V^, V h v) H i,o( Q ) 



> 



CO 



-£ft(v) 



L2(Q) + ^ 



2 



L 2 (Q) 



d-1 

E 



^2 (cV, A n ( V^Mft (t ) ) VftU, V/A . w) L 2 (n) 



> (|-^o 



> — 



2 



(3.36) 



for all w G i7 1 ' 1 (n) d if < R < R for some R G (0, 1] sufficiently small. 
To obtain higher regularity, we will use: 

Lemma 3.5 There are constants Co > 0,Ro £ (0, 1] independent of n G N U {oo} 
(and R G (0, R }) such that, if w G H^ 2 {n) d n H 3 >°(Q) d solves 

—^div h (A n (V h u h (t))V h w) = f inV'(n) 

for some f G iS" 1,0 (fi) and < /i < 1 and P n V h u satisfies (3.32) for < R < R , 
then we have 



V±e h (w),V 2 w 



< Co \\\h 2 f\\ H i,o {n) + 



\e h (w) 



H 2 >°(n) 



(3.37) 
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// additionally 
then 



e d • A n (V h u h (t))V h w\ ±i = 0, 



max 

3=0,1 



< Co||/||i?1.0(Q). 



(3.38) 



(3.39) 



Proof: Let < Rq < 1 be at least as small as above. First of all, since A n (0) = 
D 2 W(0) , we obtain 

div h (A n {0)V h w) = ^d Xd (D 2 W(0)V h w) d + dw xl (D 2 W(0)V h w)' 
= ^(D 2 W(0)8 2 Xd w e d ) d + ^(D 2 W(0)d Xd (V x ,, 0)w) d + dw x ,(D 2 W(0)V h w)' 

where A' = {aij)i = \^„ d ,j=i,...d-i for A G lR (ixd . We note that the second and third 
term consists of terms of V x /VhW. Moreover, 



{D 2 W{0)d 2 w®e d ) d = Mdl d w 



for some symmetric positive definite matrix M , which follows from the Legendre- 
Hadamard condition (3.5). Hence 

^d 2 Xd w = M- 1 (div ft (QV h w) - ^(Qd Xd (V x ,,0)w) d - dhv(QV^)') 
for Q = D 2 W(0) and therefore 



1 



< C 



div h (D 2 W(0)V h w 



+ \\V x >Vhw\\ H i, 



_ff!.0(Q) 



Thus Korn's inequality and \\d Xd ^e h (w)\\ L 2^ < \\Vlw\\ L 2^ yields 



v\e h {w),V 2 h w 



< C 



H" 1 '°(f2) 

d\w h (p 2 W{V)V h w) 



+ 



(3.40) 



Next we use that 
di\ h {A n (V h Uh)V hw) 



= di Vfe (D 2 w{Q)v hW ) + J d\v h (p n D 3 w(TPZv h u h )ip!;v h u h ,p!;v h w}) d T 

= div h (d 2 W{0)V h w) + d\v h (G(PZV h u h )[P%V h u h , P^ h w\) , 
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where G G C°°(5 £ (0); C 3 (R dxd )) for some suitable e > 0. Hence Corollary 2.5 
implies 

\\G(P%V h u h )[PZv h u h ,PZV h w]\\v h 

< C\\G(P%V hUh )\\ Vh \\PZV h u h \\ Vh \\PZV h w\\ Vh 

< Ci2o||(V ft ti;,V^)|| Hl , . 

where \\f\\v h = ||(/) Vfc/)||#i,o and we have used (3.11) as well as (3.32). Hence 

4 iVh (n*w(°)v^)\\ HV(a) 

< \\dw h (A n {Vu h )V h w) || H i,o ( n) 

V h (G(P%V h u h )[Py h u h , P%V h w] 



+ 



^ \\ h2 f\\m<o{a) +CIt * 



e h (w),V 2 h w 



(3.41) 



Combining the latter estimate with (3.40) for sufficiently small Rq £ (0, 1] , we obtain 
(3.37). 

Now, if additionally (3.38), then 



j^(A n (V h u h )V h w,V h p) H i,o = (f,<p) H i,o 

for all ip e V(Q) d . Hence, choosing tp = d 2 Jwo with wq = w — p| f Q w dx and 
I7I < 1 and using integration by parts, we obtain by (3.36), (3.37), and (3.42) below 



sup 



2 



< C ||/||Hi,o ( n)max 

< Co[||/|| H i,o (n) +i2 



d 2 7 



w 



+ CR 



h 



£h(w) 



h £h{w) 



J max 



1 



i? 1 '°(n) 



Thus, choosing Rq sufficiently small, we obtain 



< Co || /|| ifi.o^) 



with Co > depending only on Q . This finishes the proof. 



Lemma 3.6 Let P n Vu h (t) satisfy (3.31) for some n <E N U {00}, < h < 1, 
t 6 [0, T] , and < R < Ro , where Ro € (0, 1] is so small that all previous conditions 
are satisfied. Then 



(id^A n (V h u h (t)))V h w,V hV y 



< CR 



£h(w) 



£h(v) 



L 2 (Q) 



(3.42) 
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if 1 < \P\ < 2 and 



((afA„(V^(t)))V^,V^) 
i/ \/3\ = 3 and (3^ 3e , i.e., fl£ / <9 t 3 . Moreover, 



















^) 






V(fi) 


L 2 (n) 



(3.43) 



^2 ((dA^(V^(t)))V ft w,V^) n 



< CR 



h £hH 



- h e h {v) 



L 2 (n) 



(3.44) 



for all j = 1, . . . , d — 1 . T/ie constants C are independent of V/jU/j(t), w, v, h,n, R. 

Proof: If \p\ = 1, then (3.42) is just (3.35). Next let \f3\ = 2. Then for j,k = 
0,...,d-l 

0*A fc 4,(V ft u h ) = Pn J D 3 W ? (PnV^)(P n a Zj 5 Zfc V fe n /l ) 

+P n Z) 4 W ? (P n V ft n /l )(P n a..V^,P n 5 Zfe V h n /l ), 



where 



P n d Zj d Zk -e h (u h ) 



< C 



H 1 .o(n) 



P n V z -e h (u h ) 



< C Rh 



V(«) 



due to (3.31). Together with (2.15) the latter estimate implies (3.42) in the case 
\(3\ = 2. Moreover, (3.44) is proved in the same way using that dtd Xj \eh(v>h) G 
H 1 ^) is uniformly bounded and again (2.15). 
Finally, if \(3\ = 3 with d£ / d t 3 , we use that 

d Zj d Zk d Zl A n (V h u h ) = P n D 3 W(P n V h u h )[P n d Zj d Zk d Zl V h u h ] 
+P n D 4 W(P n V h u h )[P n d Zj 8 Zl V h u h , P n d Zk V h u h ] 

+p n D 4 w(p n v h u h )[p n d Zj v h u h ,p n d Zk d Zl v hUh ] 
+p n D 4 w(p n v h u h )[p n d Zt v h u h ,p n d Zj d Zk v hUh ] 

+P n D 5 W(P n V h u h )[P n d Zj V h u h , P n d Zk V h u h , P n d Zl V h u h ] 

Since P n V z V h u h G L°°(ft) and P n V 2 z V h u h G Jf 1 - ^) L 4 ' 2 (ft) are of order 
due to (3.31), the estimates of all parts in 



^ ^A n (V h u h ))V h w,V hV y 



h 2 

which come from terms involving D^W or D 5 W can be done in a straight forward 
manner by 



CR 



\e h {w) 



L 4,oo 



< C'R 



L 2 



\e h {w) 



L 2 
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uniformly in < h < 1 and n £ No U {oo}. It only remains to estimate the part 
involving the D 3 W -term: To this end we use that (3.31) and (2.14) imply 



^ ((D 3 W(P n V h u h ))[P n dPv h u h , P n V h w],P n V h v} 



< 9* 
- h 



< CR 



n 



h 



V(ft) 



L 2 (U) 
1 
h 



-e h {w),V h w 



V(Q) 



L 2 (Q) 



L 2 {n) 



Altogether we obtain (3.43). 



Corollary 3.7 Let P n V h u h (t), n G N U {oo}, < h < I, < R < R be as in 
Lemma 3.6. Then we have 



1 



^2 ((d t A n (V h u h (t)))V h w, V h v) H i, 



1 



^2 ((d?A n (V h u h ))V h w,V h v) Hlfi 



< CR 


\e h (w) 




\e h {v) 


< CR 


\e h (w) 


V(Q) 


\e h {v) 



H 1 ' 



uniformly in < h < 1 , < R< R , ra G N U {oo} , t £ [0, T] , and < T < oo . 
Proof: Because of (3.35), we have 



1 



((dtA n (y h u h ))V h Wj, VftVj-) £2( n) 



< CR 



h 



Sh{Wj) 



L 2 



e h (vj) 



L 2 



for j = 0, . . . ,d-l and (w ,v ) = (w,v), (vjj,Vj) = (d Xj w,d Xj v) if j = 1, . . . , d- 1 . 
Moreover, 



1 



-2 ((^^^(VftUh^Vftt^VfcUj) 



/? 2 



L 2 (n) 



h 



£h(w) 



h £h{vj ' 



< C'R 



L 2 



h 



H i.o 



/l 



if 1 ." 



for j = 1, . . . , d — 1 due to (3.44). Altogether this implies the first estimate. 
Similarly, (3.42) yields 



\ {{ d t A n^hUh)WhWj,V h Vj) L 2 



/l 2 

CR 



(fi) 



< C"i2 



L 2 



H 1 ' 



where (wj,Vj) are as above. Finally, (3.43) implies 



h 



v(n) 



£h(vj) 



< C'R 



L 2 



h 



v(n) 



h 



£h(v) 
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This shows the second estimates. 



Next we will show solvability of (3.19) and the estimate (3.22) for k = 0. 

Proposition 3.8 Let < h < I, P% , n G N, be the smoothing operators from 
above, let = I, and let F n be defined as in (3.14)- Then there are constants 
C > 0,M G (0, 1] independent of n G N U {oo} such that for any f G H 2 '°(£l) d 
with \\f\\ H 2,o < M h and f n fdx = there is a solution w G H 2 ' 2 (Q) d n H A ' (Q) d , 
which is unique up to a constant, such that 



1 



for all (p G H± er {Q) d and 



^e h (w),V^e h (w),V 2 h w 



for some Co > independent of h, f,n. 



jf 2 .o(n) 



< Co||/||ij2,0(Q). 



(3.45) 



(3.46) 



Proof: 

Because of (3.14), (3.45) is equivalent to solve 
{L h w,ip) H -i jH i = (D 2 W(0)V h w,V h ^ 



L 2 (Q) 



= (/,v)L2(n) - p ( i? n(V/ l u'),V /l ^) L2OT = (G h (w), V¥>) £2( n). 



(Q) 



We will prove the proposition with the aid of the contraction mapping principle. To 
this end we note that for every / G Jf 1 >°(n) d and F G H 1 ' 1 ^** n H 2 >°(n) dxd 
there is a unique w G H 2 ' 1 (n) d n H 3 >°(n) d such that 



(L h w,ip) H -i Hl = (/,<£) L 2(Q) + (F,V^) L 2 (C) 



(3.47) 



for all 93 G Hl(Q) because of the Lemma of Lax-Milgram, Korn's inequality, and 
since L/j commutes with tangential derivatives. The solution satisfies 



< Co (||/||Hi,o ( n) + ll^ll^.o 



H 2 .o(Q) 



(O) 



(3.48) 



for some universal Co > 0. Moreover, (3.47) implies 

—Ldw h (D 2 W(0)V h w) =f- dw h F in V'(Q). 

Therefore w G H 2,1 (Q) d by standard elliptic regularity. Hence Lemma 3.5 together 
with the previous estimate imply 



-e h {w),V-£ h {w),V 2 h w 



< (||(/^ 2 V fc F)|| Hl ,o (n) + 11^11^.0^) 
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for some universal Cq > 0. Using (3.33) and Corollary 2.7, one derives that 



\\G h ( Wl )-G h (w 2 )\\ Hl ,o m <CM 



v(n) 



for some C > provided that 



max 

J"=l,2 



< 2C M h, 



(3.49) 



if 1 '°(Q) 

where Co > is as (3.48) and Mo G (0, 1] . Here we note that 

d Xk G h { Wj ) = -^A' n (V hWj )V h d Xk Wj, 
d Xk d Xl G h { Wj ) = -^A' n {^hWj)V h d Xk d XlWj 

- ^P^WiPnVhW^PnVhd^Wj^nVhd^Wj] 

for all k,l = 1,. . . ,d- 1, j = 1,2, where ^(VfciOj) = A n (V h Wj) - D 2 W(0). To 
estimate the A' n -terms one uses (2.14) or (3.33) and to estimate the D 3 Ty-term one 
uses (2.15). 

Furthermore, using Corollary 2.5, one shows in the same way as in the proof of 
Lemma 3.5, that 

h 2 1| V h (G h ( Wl ) -G h (w 2 )) || H i,o (a) < CM \\(Vl( Wl - w 2 ),V h ( Wl - w 2 ))\\ Hlfi 

for some C > provided that (3.49) holds. Hence, if M$ G (0, 1] is sufficiently 
small, we obtain that L^Gh- Xh — ► -X^ restricted to B 2 c o M o h(0) is a contraction, 
where is normed by 



w 



-e h (w),V-e h (w),V 2 h w 



i? 1 .o(n) 



Therefore we obtain a unique solution w solving (3.45) and satisfying (3.46) with 
H 2 '°(Q) replaced by if 1 '°(0). In order to obtain (3.46), one can simply use that 
Wj := d Xj w , j = 1, . . . , d — 1 , solves 

(A n (Vhw)V h Wj,Vh(p) L 2 {n) = (d Xj f,ip) L 2 (n) for all tp £ H] er (VL) 

and apply Lemma 3.5. ■ 

The following lemma contains the essential estimate for the linearized system to 
(3.1)-(3.4): 

Lemma 3.9 Let < T < oo , < h < 1, < R < R , n G N U {oo} be 
given, and let Rq be as in Lemma 3.5. Assume that Uh satisfies (3.31) and that 
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/ G Wl(0,T;H 1 >°) d , w G ^(^ntf 3 '°(ft) d , w x G F(ft) d . T/ien i/tere is a unigue 
to G C°([0,r|; fl 2 - 1 ^) n 5 3 >°(n)) d n ^([O,!]; ff 1 - (fi)) <i tfurf so/ves 



- ^div ft (^ n (V/ l «/ l )V/ l u; 



-4n(V A u /i )V ?i ti; e d \ Xd=± i 

M Xj =L 

{w,d t w)\ t=0 



f 





-L > 



(«j , u>i). 



(3.50) 
(3.51) 
(3.52) 
(3.53) 



Moreover, there are some constants Cl,C > 1 depending only on 17 and W such 
that 



1 



1 



d t w, T eh{w),V x ,t-r^h{w 



h 



h 



C([0,T];H^) 



(3.54) 



< C L e 



CRT 



w[(o,T-,mfi) 



^ + 



-e h (w 1 ),w 2 ,f\t=o 



iji.o 



where 



w 2 



1 



div/ l (^„(Vhii/ l | t =o)V/ i u;o) + /|t=o- 



(3.55) 



Proof: Existence of a solution w G C 1 ([0, T]; H 1 ' (Q)) n C°([0, T]; F(fi)) can be 
obtained by the standard energy method, cf. e.g. [19, Theorem 10.8] with H = 
H 1 '\n) d ,V = V(n) d and 



a(i; v, w) 



1 



(n), 



v,w g v(ny 



which is bounded and coercive on V(f2) d because of (3.36). If n 7^ 00, then u; G 
C 2 {[0,T};H lfi (Q)) n C 1 ([0,T];V r (O)) can be obtained by the same technique as in 
[22, Section 30.1], where we note that 

^div h ((d t A n (V h u h ))V h -): V(n) - H^°(Sl) 

is a bounded linear operator with operator norm bounded uniformly in t G [0, T] 
because of the smoothing operator P n in the definition of A n and dt^h u h G 
C([0,T];L 2 (ft)). Moreover, w G C°([0, T] ; i/ 1 ' 2 (fi) n H 3 '°(Q))) follows from (3.50) 
with dfwj G C ([0,T];F 1 '°(^)), (3.28) with n < 00 and standard elliptic the- 
ory. Finally, if n = 00, then existence of a solution u; G C 2 ([0, T\\ if 1,0 (f2)) n 
C°([0,T]; J ff 1 ' 2 (fi) n i7 3 '°(17)) can be obtained from the case n G N by using the 
uniform bounds due to (3.54) proved below and passing to the limit n — > 00. 

Hence the main task is to establish (3.54). First of all, we note that (3.50)-(3.52) 
imply 

a(t) := / w(t)dx= / w$dx + t I w\dx+ / (t — r) / f(T,x)dxdr. 
Jn Jn Jn Jo Jn 
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Hence, replacing w(t) by w(t) — a(t) and subtracting from (wo,w\,f) their mean 
values with respect to , we can reduce to the case 



wq dx = / w\ dx 



n 



/ f(t)dx = / w(t)dx = 



for all < t < T. 

Now we differentiate (3.50) with respect to t and multiply with d^w in H := 
H 1 ' ^). Then we obtain 



1 d 



~ ( ||9 t 2 w||^ + ^ {A n (V h u h )V h d t w,V h d t w) H 



1 



< |(d t /,3 t yy + 



((d t A n (y h u h ))V h d t w,V h d t w) H 



^2 ((^ 2 ^„(V/ l u fe ))V/ l u;, V h d t w) H - ({d t A n (V h u h ))V h w,V h d t w) H 
in the sense of distributions, where we have used 



_d / 1 



1 



dt 



\2h? ( A ^( V hUh)^hd t w,V h d t w) H + p {{d t A n {V h u h ))V h w,V h d t w) H 



1 d 



1 



+p ((9 2 ^„(V^))V^,V^) H 
and (3.51)-(3.52). Due to Corollary 3.7 we have 

-^\((d t A n (V h u h ))V h d t w,V h d t w) H \ < CR 

-^\((d?A n (V h u h ))V h w,V h d t w) H \ < CR 



(3.56) 



-e h (d t w) 



h 



V(Q) 



£h{dtw) 



ii 



for every t € [0, T] . Moreover, because of Corollary 3.7 again, 



sup 

0<T<t 



j^((d t A n (V h u h {T)))V h w(T),V h d t w{T)) H 



< CR 



h 



£h{w) 



L°°(0,t;H) 



£h(d t w) 



L°°(0,t;H) 



Therefore the previous estimates, (3.36), and Young's inequality imply 

2 



sup 

0<T<t 



1 



d t w(r), T e h (d t w(T)) 



h 



< CR 



+C 



d^w, -^e h (w),V XttT e h (w 



e h (w 1 ),w 2 



h 



+ CR 



H 



L 2 (0,t;H) 
2 



c o\\dtf\\h(o,T;H) 
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Now 



h 



£h{w) 



L°°(0,t;H) 



< Cn max 
Vi=o,i 



-e h (d{w) 



+ 



L 2 (0,t;H) 



-£h(wo) 



H , 



due to 



L°°(0,t;H) < C (\\f\\ W i(0,t;H) + ll/|t=o||ff) 



(3.57) 



with some Co > independent of t > 0, cf. (2.4), and 



^e h (w),V^£ h (w] 



L°°(0,t;H) 
i2„ 



< Co (\\fh°°(0,t;H) + \\dt w \\L^(0,t;H)) 

< C (\\f\\wl(0,t;H) + ll/|t=0||ff + \\dtW\\L^(Q,t;H)) 

due to (3.39) and (3.57) uniformly in < t < T. Hence we conclude 



sup 

0<T<t 



< CR 



9fW, ^e h (w), V Xjt ^£h{w)) 



H 



L 2 (0,t;H) 



+Co\\f\\wi ( o,T;H) +C 



£h(wi),W2,f\t=o 



where we have used R < 1 and (3.55). Therefore the Lemma of Gronwall yields 

2 



L°e(<d,T;H) 



< C L e 



This shows (3.54). 



CRT 



h 



£h(wi),w 2 J\t=o 



+ 



H 



WfWwl(0,T; 



H) ■ 



Finally, we consider (3.24)-(3.27) with / replaced by — div^/ in its weak form, 
namely: 



-(d t w,d t f)Q T + -^{A n {V h u h )V h w,V h (p)Q T 

= (f,V h ip)Q T + {wi,y\ t=Q ) H -l H l 

h ' h 



w 



Xj=L 



w\ Xj =-L for all j = 1, . . . , d — 1 
w\ t =o = w - 

for all <p G C 1 ([0,T];fl^ r (fi) d ) with ip\ t=T = 0. 
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(3.58) 
(3.59) 
(3.60) 



Lemma 3.10 Assume that satisfies (3.31) with R G (0, Rq] and some given 
< h < I, n G NU{cx)} and /ei R G (0, 1] 6e so Mail i/iai (3.35) and (3.34) /io/d. 
Lei io G H l (Q T ) d be a solution of (3.58)-(3.60) for some f G L l (0, T; L?(VL) dxd ) , 
wq G L 2 (Q) d , and w\ G Hp er (Q) d and let u(t) = f*w(T)dT. Then there are some 
Cq, C > independent of w and < T < oo snc/i t/iai 



(3.61) 



11 / L^(0,T;L 2 ) 

< C e CRT (||/|| L l (0 ,T;L2) + IK|| L 2 + \\wi\\ H -i (Q) ) ■ 



rpf 

Proof: Let < T < T and define v,T'(t) = — L w(r)dT. We choose <p = 
ut'X[o,t'] m (3.58) (after a standard approximation). Then 

\\\w{T')\\ 2 LHn) + ^(A n (V h u h )V h u T ,(0), V h «rK0)) n 
= --^((dtA n (V h u h ))V h u T ',V h u T ')Q T , - (/, V h tt T ')Q T / 
-<wi,«r'(0)) H -i |H i + ^INHi2 (Q) 
Hence (3.34), (3.35), and u T ,(0) = -u{T') imply 



HT')\\ 2 L2{n) + 



< CR 



L 2 



+ c \\ w o\\l 2 (n) + c (|I/IIli(o,t ; l2) + IHUff- 1 ) 



L 2 (Q T i) 
1 
h 



£h(u) 



L°°(0,T';L 2 ) 



for all < T' < T. Since u T -(t) = -u(T') + u(t), we obtain 



£h(uT') 



< 



L 2 (Q Tl ) 



h 



+ T' 



L 2 {Q T i) 



h 



e h (u(T')) 



L 2 (n) 



Hence there is some k > independent of R G (0, Rq] , h G (0, 1] , such that 

2 



MIl°°(0,T';Z, 2 ) 



+ 



< CR 



h £h{u) 



\e h (u) 



L 2 (Q Tl ) 



L°°(0,T';L 2 ) 

+ Co (||wo||£ 2( n) + Ikillff-^n) + H/ll£i(o,T;L*)) 



provided that RT' < n. By the lemma of Gronwall we obtain (3.61) for all < 
T < oo such that RT < k. Now, if < T < oo with RT > k, we apply the 
latter estimate successively for some = To < T\ < ... < T/v = T such that 
R(T j+1 - Tj) < k, j = 0, . . . , TV - 1, and JV < 2Rk~ 1 T . Hence we obtain 

w, \e h {u) 

n / L°°(0,T;L 2 ) 

< (C ) N e CRT (||/|| L1(0 ,T;L 2 ) + IKH L 2 + lki|lH-i ( n)) > 
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where 

(Co)* < exp (2k" 1 InCoRT) < exp(C i?T) 

since N < 2Rk~ 1 T . This implies (3.61) for some modified Cq,C independent of 
R G (0,.Ro], h G (0,1], < T < oo. ■ 



3.4 Local in Time Existence 

For the following we assume that 9 > 0, < T < 1, and («2,h, "3,/i)> fh are as in 
Theorem 3.3 and set w\ = u^hi w 2 = U3,h- Moreover, we assume that Rq G (0, 1] is 
so small that all the statements in Section 3.3 are applicable. - Note that T < 1 is 
not a restriction for the proof of Theorem 3.3 and Theorem 3.1. By a simple scaling 
with T -1 in time t and h we can always reduce to this case changing M > 
by a certain factor depending on T if necessary. (Of course this finally influence 
the smallness assumption of ho > in the case 9 > and the starting smallness 
assumption on M if 9 = 0.) 

Moreover, let Cl > 1 be the constant in Lemma 3.9 and let Co > 1 be as in 
(3.22). Then (3.6)-(3.7), (3.22) imply 



h 1+e d t f h 



+ max 

M<i 



+ max 
wl{[o,T]-,m.o) fc=o,i 



.1 



^h(ul h ),V T e h (ul h ),Viul h 



h 



-e h (w 1 ),w 2 , h 1+e d]f h \ t =o 



< Mh 1+e 



(3.62) 



where M = (2 + Co)M . If 9 > 0, we can find some ho G (0, 1] (depending on M) 
such that R := 6C L Mh e < R . If 6 = 0, we assume that M > is so small that 
R := 6ClM < Ro . In this case we set ho = I. 
Under the latter assumptions we will prove: 



Theorem 3.11 For every n G N and < h < ho there is some To > and a 
unique u n h G C 3 ([0, T'];H lfi ) flC([0, T];H 2 > 2 OH 4 ' ) solving (3.15)-(3.18) on (0, T') 
with T' = min(T, To) , satisfying 



max 



||(^«,i^(5X),V^^^<))|| c([0T , ];H1 , 0) < 2,C L Mh l - 



(3.63) 



uniformly in < h < ho . Here To > depends only on M, n . 

We solve (3.15)-(3.18) by solving (3.24)-(3.27). To this end, let id 
w!H +1 , k G No , be defined recursively by 



° = and let 



« +1 



^div fe (An(V ' h u k n )V ' h w k n +l 
A n (V h u k n )V h w k n +1 e d 

wt +1 



Xd=±h 



(w k+ \d t w k+i : 



t=0 



d t f h h 1+e 
0, 

w k+1 

Xj =L 
(W0,n,wi), 



in n x (0,T'), 



,d 



(3.64) 
(3.65) 

(3.66) 
(3.67) 
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where w 0>n = u™ h , 

u k (x, t) = vJ Q \{x) + f w k (x, t) dr for all k G N , (3.68) 
J o 

and v>Q h ,u™ h solve (3. 19), (3. 20), respectively, cf. Remark 3.2. The existence of a 
unique solution wf^ 1 follows from Lemma 3.9. 

As usual for short time existence of hyperbolic equations, we first show bound- 
edness of (u^keN in some suitable "high norms" and then convergence of (ufykeN 
in some "low norms" provided that T = T (n) > is sufficiently small and n£N 
is fixed. 

In order to get "Boundedness in High Norms", we show that (w k )k£N satisfies for 
sufficiently small T' = T'(n) G (0, T] 



d 2 w k n ,-e h (w k ),V x , t -e h {w k n , 



< 2C L Mh 1+e 



(3.69) 



uniformly in < h < ho, where Cl is the constant in Lemma 3.9. To this end we 
use: 

Lemma 3.12 There is some < T'(n) < min(l,T) depending only on M,n and 
min(l,T) such that, if u^w^ satisfy (3.68) and 



dfw k ,^e h (w k ),V x ^e h {w k ) 



< 4C L Mh 1+e 



C([0,T>];Hi.°) 

then the solution w k+l of (3.64)-(3.67) satisfies 



(3.70) 



d^w k+ \le h (w k+1 ),V x ^e h (w k n +1 ) 



< 2C L Mh 1+e 



c([o,T'] ; *ri.°) 



and u k satisfies 



max 

ItI<i 



1 



+ 



d2P n V h u k ,d]P n -e h (u k ] 



d?p n v h u k ,d?p n y h (u k ) 



L°°(0,T';V) 



L^^T'-H 1 ' ) 



(3.71) 



< 6Ci C L Mh l+e <CiR h 



uniformly in < h < ho , 



Proof: First of all, if (3.70) holds for some k G N and some T' > 0, then P n VhU k 
satisfies 



d t P n V h u k ,d t P n -e h (u k ) 



C([0,T'];V) 



+ 



1 



d 2 t P n V h u k n ,d 2 t P n -e h (u k n , 



< 4C x C L Mh x + 6 



C([0,T'];H1.0) 
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uniformly in < h < ho where C\ > 1 can be chosen as the same constant as 
in (3.31). Moreover, there is some C n > depending only on n G N such that 
\\Pnf\\m(n) < C n ||/|| L2( n) uniformly in < h < 1, cf. (3.12). Now, using 



V h u k n = V h ul h + j\ h w k n {T)dr, 
Jo 

(3.22), (3.70), and the previous estimate, we conclude 



(3.72) 



max 

|7|<1 



dlP n V h u k n ,d2P n -e h (u k n ] 



+ 



diP n V h utd 2 t P n ^e h (u" n 



C([0,T'];V) 
C([0,T'];ffi.O) 



(3.73) 

< (4dC L + 1 + C' n T')Mh 1+e 



uniformly in < h < ho, < T 1 < T , j = 1, . . . , d — 1 for some C' n depending on 
n G N. Hence, if V = T'(n) > is so small that 1 + C' n T < 2dC L , we obtain 
(3.71) uniformly in < h < h , where QC L Mhl < R by the choice of h ,M . 

Because of (3.71), P n V h u k satisfies (3.31) for T replaced by V = T'(n) and 
given n G N. Hence we can apply Lemma 3.9 to conclude 



^ r 6C'Mhf.T> ( ,1+6»|| f || 

< c L e \\h\\w?(0,T>;m,°) 



CQ0,T'],m.°) 



+ 



-e h (w 1 ),w 2 ,h 1+e d t f h \t=o 



< C L e c ' RoT 'Mh 1+e 



(3.74) 



for all V > such that l+C' n T' < 2C X C L . Hence, if T'(n) > is chosen sufficiently 
small, we obtain C^e R ° T < 2Cl ■ This shows the estimate for w k+1 . ■ 

Proof of Theorem 3.11: Let , be defined as above. Because of the latter 
lemma (3.69) holds for all k G N provided that T' = T'{n) > is chosen as in the 
lemma. 

In order to show "contraction in low norms", let z k+1 : = w k+l — w k and Z k+1 : = 
u k t +1 — u^. Because of (3.64)-(3.65), z k+1 solves 



-{d t z k+ \d t ip) QTl + -.(^(VftOv^W^o 



l 



fc+i 



h? 



T' 



= -tf(MV h u k ) - A n {V h u k n - l ))V h w k n , V h <p)Q T , 
for all <p G C\[0X]lL 2 (n)) d nCP([0,T , ]-,H 1 (Sl)) d with ip\ t=T = and (z k+1 , d t z k+l ) 



t=o 
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(0,0). Here 



1 



^((A»(V h <) - A n (V h u k n 1 ))V ft <, V h tp)Q T , 



< 



CT' 



P n V h w k ,P n -e h {w k ) 













) 












L°°(0,T';V) 


L°°(0,T';L 2 ) 



L°°(0,T 



< C'T'Mhl 



L°°(0,T';L 2 ) 

because of Corollary 2.7, (3.70), and 



h £M 



L°°(0,T';L 2 ) 



1 



{(A n (V h u«) - AniVhU^VhW^Vwh 



D 6 W{{1 - T)V h u«- L + TV h u K n )[P n V h Z«,P n V h w K n ],P n V hV ) Q dr. 



Hence 



< Ce c '^ h o T 'T' 



\e h {Z k ] 



C([0,T'];L 2 ) 



C([0,T'];L 2 ) 

due to Lemma 3.10. Therefore there is some Tq > depending only on M such 
that 



dtZ k+ \\e h (Z k ^ 



< 



1 



dtZ k , -e h (Z k ] 



L°°(0,T';L 2 ) 



L°°(0,T";L 2 ) 

provided that < T"(n) < Tq. This shows that (wjDfceN converges to some 
u n e C([0,T'};H l (n)) as k -> oo and d t u k -► d t « n in C([0, T']; L 2 (S7)) for some 
sufficiently small T'(n). Because of (3.69) and (3.71), ii n and := dtu n satisfy 
the same estimates as u^, w%, respectively. By interpolation V^u^ ^fc-^oo ^hU n 
in L°°(Qt') strongly. Therefore u n and w n solve (3.24)-(3.27). Consequently u n 
solves (3.15)-(3.18). 

Finally, we have to prove (3.63). We know that u n constructed above on 
(0,T'(n)) satisfies 



, 1 1 

dfu n , d t -e h (u n ),d t V x ,tTEh(ur, 



C([0,r'];/fi.O) 



h " v lln 1 

due to (3.69). Moreover, using (3.72) and (3.62), we conclude 
1 , , 1 



< 2C L Mh l+e 



max 

3=0,1 



< 3C L Mh 1+e . (3.75) 

C([o,T']-m,o) 

Hence it remains to show that we can replace one time derivative in (3.75) by d Xj , 
j = 1, . . . , d — 1 . To this end let w J n = d Xj u n , j = 1, . . . , d — 1 . Then wP n solves 



d 2 t w{ - -^dw h A n (V h u n )V h wi = d Xj f h h 1+e in Q x (0, T') 



A n {V h u n )V h w 3 n e d \ x 



wi 



0, 



wi 



, j = 1, . . . ,d - 1, 



«,dtwi)\ t=Q = K )h X,ft) 



(3.76) 
(3.77) 
(3.78) 
(3.79) 
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with u? kh = d Xj u% h , k = 0,1. Hence Lemma 3.12 with 



w 



k _ 



k _ 



W 



k+l 



w 3 n, and dtf replaced by d Xj f implies the estimates of the remaining terms 



in (3.63). 



3.5 Uniform bounds and Proof of Theorem 3.1 

We know that u n constructed above on (0, T'(n)) satisfies 



max 

| 7 |<1 



1 



1 



d t d]u n , d2-re h (u n ), d]V Xt t-r£h{u n 



h 



h 



C([0,T'];H1.0) 



< 3C L Mh 1+e (3.80) 



where M = CqM + 2M due to Theorem 3.11. In the case 9 > 0, it remains to 
show that we can replace T'(n) by an arbitrary < T < 1 if < h < ho for 
some sufficiently small < ho < 1 independent of n G N. If 9 = 0, we will show 
that T'{n) can be replaced by < T < 1 if M is sufficiently small (independent of 
n 6 No). 

To this end we apply Lemma 3.9 for w J n = d Xj u n and w n = dtu n using the 

equations for wii,w n and that (3.30) is valid for u = u n and R = 6CiMh d < Ro, 
which implies (3.31) uniformly in n G No U {oo}. Thus we obtain 



max 

|7|=1 



d?d]u n , d]-£ h (u n ), V x , t d2-£ h (u 



h 



h 



< C L e c '^ e Mh 1+e 



uniformly in n G N and < h< ho. Due to (3.22) and (3.72), we conclude 



max 



1 



1 



d t d2u n , d]-e h (u n ),V Xtt d]-e h (u n ) 



h " v ^ *h 
If 6 > 0, we can now choose < ho < 1 so small that 

l2 



< 2C L e c '^ e Mh 1 + e 



C([0,T'];/fi.0) 



max 

It|<i 



<)f<nu n ,d]^e h (u n ), V x , t d]^e h (u r . 



< 3C L Mh 1+e 



C([0,T'];H1.0) 



uniformly in n G N and < h < h . If 9 = 0, then we choose M = (2 + C )M 
sufficiently small to obtain the same estimates. Hence we get a uniform bound as long 
as the solutions exists. Applying the results of Section 3.4 with the only difference 
that the bound M of the norm of the initial values is replaced by 3ClM and a shift in 
time by T'{n) , we can continue the solution Uh on [0, T] , T = min (T'(n) + T" , T) , 
where T" > depends only on n and 3ClM . This extended solution can be 
constructed such that (3.70) is valid. But the arguments of this section show that 
then even 



max 



1 



1 



d t d]u n , d]-E h (u n ),X7 X)t d] T e h (u, 



h 



< 3C L Mh 1+e . 



C([0,T];tfi.°) 



Therefore we can repeat this argument finitely many times and can continue the 
solution for all t G [0, T] . 

Using the uniform bounds, it is easy to pass to the limit n — > oo for a suitable 
subsequence. Uniqueness of the solution can be shown by similar estimates as in the 
proof of Theorem 3.11. This finishes the proof of Theorem 3.1. 
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4 First Order Asymptotics 

Throughout this section we assume that 

for some given g G W 2 (0, T; H* er ((-L, L)"" 1 )) n ^(0, T; H* er ((-L, L)"" 1 )) . As 
seen in the proof of Lemma 3.5, we can assume without loss of generality that 
/(_£ g{x')dx' = 0. Moreover, we assume that < < 1. 

For simplicity let = dist(F, SO(d)) 2 , which implies D 2 W(0).F = symF. 

Let L2V = divV(v) and let v be the solution of the (d — 1) -dimensional wave 
equation 

d 2 v + ^A 2 x ,v = g in (-L,^ 1 x (0,T), 



v\ Xj =-L = v\ Xj=L for j = 1, . . . , d - 1, 
(u,e%u)|t=o = 0o,^i) in 



where t> G H£ er ((-L, L) d ~ 1 ),v 1 G Hp er ((-L, L) d ~ l ) . By standard methods the 
latter system possesses a unique solution 

«GC 2 ([0,r];^ r ((-L,L) d - 1 ))nC 1 ([0,T];ij; er ((-L,L)' i - 1 )) 
nC7 ([0 > r];flJ !r ((-L,L) d - 1 )). 

Using w, we define an approximate solution u h of (3.1)-(3.4) by 



Then 





(l8 X d ~~ 24 X d ~~ 24lB)^x' V 



V (24 





24 X d - ^)A 2 ,1, 



(is x d ~ ~k x d ~ 2TT6) v x'^x' u 



Hence 



e h {u h ) = h 2 + e 



-x d V 2 x ,v %{x d -\)V xl A x ,v 

k( T 2 _ l\T7T A .„, /,2c 1 „ . _ 1^ A2 



+ j,4+e ( (14 ~ \x d )V 2 x ,A x 'V 

\2V48 X d 24 X rf 24-16/ V rr'^-\r'< 



2(48 X rf 24 X d 24-16^ %'^x' V 
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and therefore 



Moreover, 



\e h {u h )e d \ Xd=±h = 0. 



-^dw h e h (u h ) 



h 



-x d V x 'A x iv + XdV x >A x >v 



(4.1) 



, h 2+e ( ( l A ~ \x d )V x >A 2 x ,v + \{±x d - \xl)V x ,Al,v\ 

+n \ h(± r 2 _ J_ r 4 MA 3 V I 



= h L 





T2 Ar' 



where 



\\ r h\\c([0,T];L 2 (Q)) < Ch 2+6 . 

Thus Uh is a solution of 

d 2 u h - ^d\v h (D 2 W(0)V h u h ) = f h h 1+e - r h 



(d 2 w(o)v 



hUh e d 







Uh\ Xj=L = u h \ Xj= _ L , 
(uh,d t Uh)\ t =o = (uo,h,ui :h ), 



j = l,...,d-l, 



(4.2) 



inftx(0,T) (4.3) 



where 



u 



. h{x) = h i+e ^0 \ + h2+e f-x d V x ,v^ + ^ 4+e |^z 3 - \x d )V x ,A x ,v^ 
+h 5+e 



J_ T 2 _ J_ 4 MA 2 

48 X (2 24 x d 24-16/ ^x' U J 



3 =0,1. 



We will compare this approximate solution with an exact solution of the d- 
dimensional system (3.1)-(3.4) for an appropriate choice of initial values. 

Theorem 4.1 Let < 6 < 1, let vo,vi,f, uq^,u\^ and Uh be defined as above. 
Then for some sufficiently small ho £ (0,1] there are initial values (uq^jUi^) sat- 
isfying (3.6) and such that 



max 
J"=0,1 



-g£h{uj, h ) - ^e h (u j>h ] 



< Ch l+2e . 



Moreover, if Uh is the solution of (3.1)-(3.4) due to Theorem 3.1, then 
dt{u h -u h ),\e h {u h -u h y\ <Ch l+2e 

n / L°°(0,T;L 2 ) 

for all < h < ho and some C > independent of h . 
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Proof: We construct the initial values u h,ui h as solution of 



-^(DW(V h u ,h),V h <p)n = (u2, h , <p)n ~ (h 1+e f\ t=0 , <p) n , 



(D A W{V h u 0A )V h u lih ,V h ip) Q = (u 3A ,ip) n - (h 1+9 d t f\ t=0 ,<ph 



/i 2 



for all ip £ Hp er ({l) d , where 



U2+j,h = h 



1+0 



-±Al,div\ t=0 



+ h 2+e (-%V x/ Al,dlv\ t= o 



J = 0,1. 



-x d V 2 xl Al, Vj 




We note that f n U2+j t h dx = and 

\e h {u 2+hh ) = h l+e (~ 
by a similar calculation as above. In particular, this implies 

< Ch 1+e , j = 0, 1, 



^h{u 2 +j,h),^£h(dif\t=o) 



where we note that / is independent of x^. Because of Proposition 3.8, the Lemma 
of Lax-Milgram, and Lemma 3.5, (uo t h,ui,h) exist for all < h < ho if ho £ (0, 1] 
is sufficiently small and satisfy 



< Ch 1+e . 



H 2 >°(Sl) 



Altogether, (uo,hi u i,hi u 2,hi u 3,h) satisfy (3.6) and (3.8)-(3.9). Moreover, we note 
that 

1 , . 1 



max 
j"=o,i 



jfh{uj,h) - ^h{Uj,h) 



< Ch 1+2e 



L 2 (Q) 



since 



jp(£h(uo,h-uo,h),£h(<p))n = ^2{F'(V h u ^ h ),\7 h (p)n + (r 0jh ,<p)n, 

^(e h {u hh -u lyh ),e h (^))n = ^{(DW 2 (V h u , h ) - DW 2 {0))V h u hh , V h <p)n 

+( r i,h,<P)n 

due to (4.3), where F' is as in (3.13) and ||(ro,h>rL,fc)||L 2 («) ^ Ch 2+e . Here we have 
used that 



h 2 



{F'(V h u ,h),V h v) 



n 



h 2 



((DW 2 (y h u , h ) - DW 2 (0))V h u 1 , h ,V h i P ) i 



< Ch 1+2e 

< Ch 1+2e 



L 2 (f2) 

L 2 (n) 



for all (p e H} )er (Q) d because of (3.33) and the estimates for uo : h,ui,h- 
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Now let Uh be the solution of (3.1)-(3.4) due to Theorem 3.1 and consider Wh = 
dtUh — dtith- Then Wh solves 

-(d t w h ,d t (p)Q T + ^(DW(y h u h )V h w h ,V h (p)Q T - (wi th ,(p\t=o)n 

= ~ ((DW(V h u h ) - D 2 W(0))V h d t u h , V h <p) Qr - (r h , v)q t 

\,,-L = w \ Xj =-L> j = l,---,d-l, 
w\t=o = W h. 



w\ 



for all <p G C 1 ([0,T];H^ er (n) d ) with <p\ t=T = 0, where w jth = d] +j u h \ t=Q - 
\t=o, j = 0,1 and satisfies (4.2). Moreover, 



p ((DW(V h u h ) - D 2 W(0))V h d t u h , V^) s 



< 



c 













j^£h(uh) 




-e h (d t Uh) 




\e h [<p) 


L°°(0,T;L 2 ) 


L°°(0,T;V) 



L 2 (n) 



c Ch 1+2e 



h 



L 2 (n) 



due to (3.33). Hence Lemma 3.10 implies 



d t (u h - u h ), -e h (u h - u h ) 



< Ch 1+2e , 



L°°(0,T;L 2 ) 



which proves the theorem. 
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